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We present a qualitatively new approach to discover Higgs bosons of the MSSM at the LHC 
using jet substructure techniques applied to boosted Higgs decays. These techniques are ideally 
suited to the MSSM, since the lightest Higgs boson overwhelmingly decays to bb throughout the 
entire parameter space, while the heavier neutral Higgs bosons, if light enough to be produced in a 
cascade, also predominantly decay to bh. The Higgs production we consider arises from superpartner 
production where superpartners cascade decay into Higgs bosons. We study this mode of Higgs 
production for several superpartner hierarchies: mq,mg > g > mh + /i; mq,mg > g > 
m.h,H,A + a^nd mq,mg > > mh + with rUg ~ jj.. In these cascades, the Higgs bosons 

are boosted, with pr > 200 GeV a large fraction of the time. Since Higgses appear in cascades 
originating from squarks and/or gluinos, the cross section for events with at least one Higgs can be 
the same order as squark/gluino production. Given 10 fb~^ of 14 TeV LHC data, with ruq < 1 TeV, 
and one of the above superpartner mass hierarchies, our estimate of S/^fB of the Higgs signal is 
sufficiently high that the bb mode can become the discovery mode of the lightest Higgs boson of the 
MSSM. 



I. INTRODUCTION 

Uncovering the origin of electroweak symmetry break- 
ing is of the utmost importance for the LHC. If the 
world is supersymmetric - in the form of the minimal 
supersymmetric standard model (MSSM) - electroweak 
symmetry breaking is accomplished through a supersym- 
metrized two-Higgs-doublet model, with couplings and 
interactions set or restricted by supersymmetry. 

It is well known that imposing the proper electroweak 
symmetry breaking minimum leaves two undetermined 
parameters in the Higgs sector at tree-level: the ratio of 
the Higgs scalar vevs, tan /?, and the mass of the CP-odd 
scalar, tua- Radiative corrections dominantly affect the 
lightest Higgs mass, raising it from the ruled-out tree- 
level value mh = Mz up to about 125 GeV for stop 
masses and mixings that do not exceed 1 TeV [T]. Decay 
rates of the Higgs bosons can also be computed largely 
independently of the details of the superpartner sector 
(so long as decays into superpartners are either kine- 
matically forbidden or rarely occur). The Higgs sector 
can thus seemingly be approximately parametrized by 
m/i, mA, tan /?. 

This has reinforced the simplification that the Higgs 
scalar sector can be searched for, discovered, or ruled 
out in isolation from the remainder of the model [2]. A 
casual glance at the ATLAS or CMS TDRs Hi] demon- 
strates this canonical view, in which discovery potential 
for the Higgs sector is plotted in the mA-tan /3 plane (with 
some additional restriction on rrih larger than the LEP 
II bound). The Achilles heal of this simplification is the 
assumption that the most promising production channels 
of the Higgs bosons are largely the same ones as in the 
Standard Model (SM). 

We demonstrate there is potentially a much superior 
way to discover Higgs bosons in the MSSM - superpart- 
ner production with superpartners that cascade decay 
into Higgs bosons. Higgs bosons from superpartner cas- 



cades is not a new idea, see e.g. [5Ul2j. but our method 
for finding and identifying Higgs bosons within the super- 
symmetric event sample is qualitatively new. We exploit 
recently developed jet substructure techniques [T3] with 
modifications that we presented in Ref. [14] to isolate the 
boosted Higgs-to-&6 signal from the Standard Model and 
supersymmetric backgrounds. The existence of a large 
supersymmetric cascade-to-Higgs rate requires relatively 
mild assumptions about the superpartner mass hierarchy. 

The notion to find and study supersymmetric signals 
through the hadronic decays of gauge bosons, as well as 
the lightest Higgs boson, was pointed out in an early use 
of jet substructure in Ref. [15]. There, however, the mo- 
tivation was not to find the Higgs, but instead to recover 
the superpartner mass spectrum using a kinematical edge 
analysis. 

In our previous paper [TJ], we pointed out that a sig- 
nal of the Higgs boson itself can often be much more eas- 
ily found within new physics, since the new physics can 
have a larger production cross section and larger fraction 
of boosted Higgs bosons. But identifying the Higgs bo- 
son in processes that invariably have busier final states 
with more jets (and potentially more hard 6-jets) required 
modest improvements to the BDRS algorithm. This is 
not unlike the situation faced by Ref. |16j in proposing a 
method to extract the Higgs signal from tth production. 

The commonalities between our previous work, 
Ref. [14], and this paper are: 

1. We seek two-body decays of a Higgs boson into hh. 

2. We use the same jet substructure algorithm to ex- 
tract this Higgs signal. 

3. The Study Points in this paper are pure MSSM. 

4. We apply fairly aggressive cuts to reduce the back- 
grounds from standard model processes. 

The main differences between Ref. [141 and this paper 
are: 



1. The LSP of the MSSM-based Study Points in this 
paper is a neutrahncj^ Our previous work, instead 
considered study points with a gravitino LSP and 
a promptly decaying Higgsino NLSP. 

2. The new physics signal is large missing energy, with 
characteristically large /f;^'^. (Our previous work, 
instead considered the new physics signal to be one 
hard 7 plus missing energy.) This means that while 
the LHC will have evidence for new physics in chan- 
nels involving large missing energy, it quite unlikely 
that the new physics signals can be readily identi- 
fied with specific processes or decays (or models, 
for that matter). 

3. In this paper, we also consider the detection of H 
and A states decaying to bb, using jet substructure 
techniques. We demonstrate that for lighter niA ^ 
200 GeV, it is possible to uncover evidence both 
both h and H/A with just 10 fb~^ of data. 

The organization of the paper is as follows: In Sec. |ll] 
we explain how Higgs bosons can be produced from spe- 
cific kinds of two-body superpartner decays. The main 
emphasis is on the qualitative features of gaugino and 
Higgsino interactions, so as to present a very clear pic- 
ture of what superpartner hierarchies provide the most 
promising source of Higgs bosons, and how the large rates 
can be easily understood. In Sec. HI we consider the 
typically largest production source of heavy gauginos, 
namely, squark production and decay into gauginos. We 
clearly demarcate which superpartners decay into which 
gauginos, so that further studies can be guided by these 
basic observations. We then consider, for specific hierar- 
chies, the prospect of finding a boosted Higgs in one of 
these supersymmetric cascades. We show that in a con- 
siderable region of the supersymmetric parameter space, 
as many as one in four typical decay chains originating 
in a squark and ending in the LSP can contain a signif- 
icantly boosted Higgs boson. In Sec. IV we show that 



the supersymmetric parameter space we consider natu- 
rally satisfies the upper bounds on the thermal relic den- 
sity. Moreover, we demonstrate how simple changes in 
the gaugino mass hierarchy (lowering Mi) can result in 
matching the cosmological density, but without signifi- 
cantly affecting the Higgs boson signal. In Sec.|V]we then 
discuss the techniques and algorithm to find Higgs decay 
using jet substructure. We compare and contrast our 
methodology with what has been used before for Stan- 
dard Model production of a Standard Model Higgs boson. 
In Sec. |VI| we present a series of Study Points that demon- 
strate the effectiveness of our algorithm in finding one or 
more Higgs bosons of the MSSM. The series of plots of 



candidate resonance jet mass are the main results of this 
paper - demonstrating that the bb mode could well be 
the discovery mode of Higgs bosons at the LHC. Finally, 
in Sec. I VIII we conclude with a discussion of our results. 



II. HIGGS FROM SUPERPARTNER DECAY: 
"GOLDSTONE REGION" 

The main focus of the paper is on Higgs bosons 
that arise from the two-body decays of neutralinos and 
charginos. 



X? ^ h/H/A + x° 
xf ^ h/H/A + xf 



(1) 
(2) 



^ It could also be a "neutralino-equivalent" , where the neutralino is 
a NLSP and the gravitino is LSP, but the lifetime of the NLSP 
is sufficiently long that its decays are not observed within the 
collider detectors. 



It is instructive to review how these decays come about, 
and why the decay rate to Higgs bosons can be sizeable 
throughout the kinematically allowed parameter space. 

The centrally important gaugino-Higgs interactions are 
the kinetic terms of the Higgs supermultiplets. They lead 
to the component interaction terms (17j 

-D^HiD^Hl - iHupHu 

-V2g' YhJHuK ~ V2gW"'H.XK + [u ^ d) , 

(3) 

where Yh^ is the hypercharge of the Higgs field. The first 
term leads to ordinary gauge boson interactions with the 
Higgs scalars, 

+{u ^ d) (4) 

The second term leads to ordinary gauge boson interac- 
tions with the Higgsinos, 

g' Yn^B^H^a'^H^ + g WI^H^f^a^^H^ + {u ^ d) . (5) 

The latter are exactly the same interactions that Stan- 
dard Model quark or lepton doublets have with Standard 
Model gauge bosons. Notice that the interaction involves 
the entire SU(2) doublets Hu,Hd, and thus, all eight real 
scalars 

Re(ff°), Im(770), ff^ , Re(i/0), Im(i/0), . (6) 

In practice, as is very well known, linear combina- 
tions of the above become the physical mass eigen- 
state Higgs bosons {h, H, A, H^) as well as the Gold- 
stone bosons {w^ , z) associated with the Standard Model 
gauge bosons , Z. The critical observation is that the 
"supersymmetrized" interactions , 

- V2g'BH^H*^ - V2gW''Hut''H* + {u ^ d) , (7) 

necessarily have gauge coupling strength to all compo- 
nents of both Higgs doublets. This means that, all other 
things considered equal, decays of B/W H + h/A/H 



2 



lead to equal branching ratios into the different compo- 
nents of the doublet. A Higgs scalar, therefore, is just as 
common as a or Z in this type of cascade decay. 

In practice, the gaugino interaction eigenstates mix 
with the Higgsino interaction eigenstates through the 
same interactions, Eq. ([t]), that led to large decay rates 
of gauginos into physical Higgs bosons. In the mass 
basis, the decays in Eqs. ([l]),([2]) roughly translate into 
heavier neutralinos Xs 4 ^-nd charginos X2 decaying into 
their lighter counterparts Xi2i Higgs bosons. 

The relevant branching ratios are thus X4 3 — > h/H/A + 
(xi or X2) ^-iid X2 ^ h/H/A + xf ■ The other possible 
decay X2 ^ ^ + Xi is mostly kinematically forbidden in 
the region of our interest. 

Another type of cascade decay occurs when the Higgsi- 
nos are heavier than the winos and/or bino. This opens 
up the decay channels H h/H/A + B/W . This might 
well provide an interesting source of Higgs bosons given a 
cascade from third generation squarks to H . Our prelimi- 
nary work on this cascade suggests it takes more luminos- 
ity than 10 fb~^ (which is the main focus of this paper), 
and requires adjustments to the cut-based search strat- 
egy to optimize for a signal of third generation squarks. 

There are yet other superpartner cascade decays that 
could also lead to Higgs bosons, such as stop decay 
t2 ^ ii + /iM/iJ [m [1^. To the extent that this pro- 
cess occurs for the specific points in the MSSM parame- 
ter space we present below, it is included in our inclusive 
analysis. In practice, however, the production cross sec- 
tion of just the heavier stop t2 is small relative to the 
large number of other squarks (and gluino), while the 
branching ratio t2 — > ti/i is also accidentally small when 
mg^ — TO(^ which we take for our Study Points. This 
implies negligible numbers of lightest Higgs bosons arise 
from t2 decay. 

Finally, there is an interesting possibility of cascade 
decays into the charged Higgs H^ . If mH± > mt + mi,, 
the hadronic decay mode iJ+ — >■ tb often has a large 
fraction, and thus could be an interesting candidate for 
jet substructure techniques, utilizing top tagging [TBI [22- 
150] or other novel methods. 



III. CASCADING TO BOOSTED HIGGS 
BOSONS 

The largest rate for Higgs boson production arises 
when first or second generation squarks cascade decay 
through gauginos, which then decay into lighter Higgsi- 
nos and Higgs bosons. Generally, first and second gener- 
ation squarks decay as 



q + W 
q + B 



(8) 
(9) 
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FIG. 1. The branching ratios of heavier gaugino-Iike neu- 
tralinos and charginos into lighter Higgsino-like ones plus the 
lightest Higgs boson is shown for the following parameters: 
We take 100 GeV < Mi = M2/2 < 400 GeV for all Fig- 
ures, = 150 GeV in plots I and III and \fj,\ = 200 GeV in 
plots II. Plots I and II have heavier sleptons, mi > 800 GeV, 
so that two-body decays are kinematically forbidden. In plot 
III, we take rrij- — 500 GeV, which allows the wino to decay 
to left-handed sleptons once M2 > 500 GeV. This is why the 
branching ratios of xl, X2 decrease above Mi/|/x| > 1.7. 



SO long as the wino and bino satisfy the simple kinemat- 
ical requirement that they are lighter than the squarks. 
The left-handed squarks can also decay to the bino, but 



this rarely happens when the wino mode is kinemati- 
cally open, since the ratio of bino to wino couplings 
for the left-handed squark doublet is proportional to 
yq^g' /gf' — 0.01. Thus, to very good accuracy, first and 
second generation left-handed squarks decay through the 
wino, right-handed squarks decay through the bino. 

Given a Higgsino plus Higgs boson lighter than the 
wino and/or bino, the two-body decays into Higgs bosons 
discussed in Sec. [TT] become applicable. Since the lighter 
quarks (w, d, c, s) have Yukawa couplings far subdomi- 
nant to the gauge coupling strengths g,g', the cascade 
in which squarks decay directly into the Higgsinos es- 
sentially never occurs. This implies the large QCD- 
dominated production cross sections of squarks can lead 
to substantial numbers of Higgs bosons from the cascade 
decays with only modest mass hierarchy requirements. 
Moreover, in addition to squark pair production, squark- 
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FIG. 2. The branching ratios for decays to the hghtest Higgs 
boson as a function of Mi / ji. The MSSM parameters for each 
plot are the same as the three rows in Fig. [l] Here g_t refers 
to the sum of ul and dz, (both components of the electroweak 
doublet), while qr refers to either ur or dR. 

gluino and gluino pair production can also lead to sub- 
stantial rates of squarks through the two-body g — qq. 

We have thus clearly demarcated the superpartner cas- 
cades into Higgs bosons as having two largely indepen- 
dent sources, namely 

qL^q + W; W ^ h/H/A + H (10) 

qn^q + B; B^h/H/A + H . (11) 

This will be useful as we consider variations of gaugino 
masses and slepton masses. 

Gauginos can also decay to sleptons, and it is fairly 
easy to see what effect they have if they are lighter than 
the bino and/or wino. For simplicity, consider all three 
generations of sleptons to be degenerate. If all sleptons 
are lighter than the bino, we can estimate the branching 
ratio by just summing over three generations of right- 
handed and left-handed leptons plus one Higgs doublet. 
We get 

BRiB^hH^)^\^^^§^^^,mHm 

(for <Ah) 



where the 1/4 comes from for picking just h from 
(/i, zq, w'^ ,w~), assuming tua is large. (It is easy to gen- 
eralize for smaller niA-) If left-handed sleptons are lighter 
than the wino, we get 

(for rrij^ < M2) . 

So, the effect of all sleptons lighter than the bino is 
to very efficiently suppress the cascade decay of right- 
handed squarks to Higgs bosons from about 1/4 to about 
a percent. Conversely, the effect of left-handed sleptons 
lighter than the wino is to reduce the branching ratio 
from about 1/4 to 1/16. While this suppression is sig- 
nificant, it certainly does not eliminate this decay mode, 
and illustrates the robustness of finding a Higgs boson 
within a fairly generic superpartner cascade. 

We can study the branching ratios in more detail nu- 
merically. As we have already seen, the likelihood of 
finding a Higgs boson in a complex decay chain originat- 
ing from a squark can be approximated, to a large extent, 
by the product: 

Pgh EE Bv{q^ X ^ h + Y) « (14) 
J2 Br{q^Xa + --.)'BT{xa^h + ...) , 

where X and Y are other particles or superpartners. 

Consider now two interesting regimes for the masses of 
the Higgs bosons. The first, "large m^", and the second, 
"smaller m^" • 



A. Large jtia 

The first regime we consider is when 

niA > min(m,j, [Ma - |Mi - /i|) , (15) 

often described as the "decoupling limit". In this limit, 
all of the Higgs bosons H,A,H^ are predominantly eigen- 
states from the second Higgs doublet. These Higgs 
bosons are much heavier than the lightest Higgs as well 
as the lighter superpartners in the model. In practice, we 
take TTiA ~ 1 TeV, and thus the scalars H, A, have 
masses ~ 1 TeV, while h mixes minimally to H with 
ordinary couplings to standard model particles. 



1. Higgs in a cascade: 

In Fig. [1] we show the branching ratios of X4 3 ^ + 
Xi,2: X2 ^ h + Xi and in Fig. [2] we show Pq^h and Pg^h 
as a function of Mi/fi in the large tjia regime. 

Plot I in Figs. |l|2| is generated with all the squarks and 
sleptons set to 1 TeV. Because of using a small value of 
\fj,\ (namely, 150 GeV), winos are relatively heavier than 
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the Higgsinos and mix minimally to the rest of the gaugi- 
nos throughout. As a result, the heavier mass eigenstates 
(X4 and X2 ) mostly winos. As indicated in Plot I in 
Fig. [1] winos decay significantly to the lightest Higgs bo- 
son. In fact, for AI2 > 300 GcV, wino decay follows the 
"Goldstone region": roughly 3/4 of the time the wino 
decays into longitudinal W/Z and 1 /4 of the time it de- 
cays into the lightest Higgs boson. For a large part of 
the parameter space the mass gaps between Xs and Xi.2 
are not large enough to allow a two-body decay into the 
lightest Higgs boson. Once outside the kinematically for- 
bidden zone, however, the branching ratio of the decay 
X3 — /i + Xi 2 rises quickly with increasing Mi . In this 
region, Mi/fi is large and X3 is mostly a bino. 

The same spectrum is used to generate plot I in Fig. [2] 
Note that the right-handed squarks decay mostly to the 
bino and so Pqi^h looks almost identical to the branch- 
ing ratio of X3 decaying to Higgs boson. Similarly, the 
left-handed squarks decay mainly to the winos and Pqj^h 
follows the partial decay width of xl and xt to Higgs 
bosons. The other feature to note in this plot is that 
Pqj^h goes down for large Mi, signifying that the decays 
of squark to quark plus wino are beginning to be affected 
by kinematical suppression from the heavy wino. 

Plots II in Figs. l|2 are similar to Plots I except that 
slightly heavier Higgsinos = 200 GeV) are used. 
Larger Mi//x is needed in order to open up X3 decays 
to Higgs bosons. A curious rise is seen in X4 decays for 
small Ml. It is an artifact of decays X4 ~^ + xf 
shutting down, thereby, causing total decay width of xl 
to shrink. This feature is more prominent for negative 
fj,. Even if all parameters in the chargino mass matrices 
are held fixed, taking M2 to have the opposite sign of /i 
reduces the splitting among the mass eigenvalues. This 
results in heavier xt and prevents the two-body decay 
X4 — >■ + xf decays even for heavier xl- Once again, 
Pq^h and Pq^^h follows the partial decay width of winos 
and bino respectively. One thing to note is that even 
though there is a sharp rise in X4 and X2^ partial widths 
for small Mi, there is no such curious feature in Pq^h- In 
this limit, M2 ~ and decays of squarks to X4 and X2^ 
suffer because of rising Higgsino content in them. 

Finally, for Plot HI in Figs. T|2 all parameters are 



the same as Plots I except that the sleptons are taken 
be lighter, in this case, 500GeV. As the wino mass is 
increased above this value, the wino-like neutralino and 
charginos begin to decay into slepton modes, reducing 
the branching fraction to the lightest Higgs boson. 



(and even with 300 GeV, as shown) . The boost was found 
by generating samples of 5000 supersymmetric events at 
different values of Mi/^i using PYTHIA v6.4 [32, and 
plotting the Higgs transverse momenta. 

Both of the plots in Fig.[3]are made with /i — 150 GeV, 
tan P — 10 and all squarks with mass of 1 TeV. Sleptons 
have mass of 1 TeV in plot I and 500 GeV in plot II. The 
presence of light sleptons reduces the fraction of super- 
symmetric production that leads to a boosted Higgs bo- 
son in the cascade. This is due not only does the overall 
lower fraction of Higgs bosons appearing in the cascades 
(see plot IIIA and IIIB in Fig. 2| when heavy neutralinos 
and charginos decay to them, but also fewer of the Higgs 
bosons in the decay chain are boosted. 
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FIG. 3. The fraction (in %) of boosted Higgs bosons as a 
function of Mi/fj, with M2 = 2Mi, p = 150 GeV and tan/3 = 
10 in samples of events generated by PYTHIA. In the plots 
the red and dotted lines represent the percentages of Higgs 
bosons with pr > 200 GeV and the green dot-dashed lines 
represent the fraction of Higgs with pr > 300 GeV. In the 
left Figure the squark masses are 1 TeV, while in the right 
Figure the squark masses are 750 GeV. All other relevant soft 
supersymmetric breaking masses are kept at or above 1 TeV. 



2. Boost of a Higgs boson in a Superpartner Cascade 

In a typical cascade, a Higgs boson appears from the 
decay of a massive superpartner. The large release in rest 
mass results in a large recoil energy, i.e., Higgs bosons 
from superpartner decays are naturally boosted. This is 
demonstrated in Fig. [3j which shows that a significant 
fraction of Higgs bosons are boosted with pT > 200 GeV 



B. Smaller niA 

The second interesting regime of the Higgs sector that 
we consider is smaller itia, where 



niA < min(|M2 - fi\, \Mi - ^| 



(16) 



There are really two distinct regimes of smaller ttia- 
the first is when all the Higgs mass eigenstates 
{h,H,A and H^) are comparable in mass and the CP 
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even neutral Higgs scalars h and H mix maximally among 
each other. The second is when there is less mixing, but 
H, A, are light enough to be kinematically accessible 
to gaugino decay. We will examine both of these cases 
below. 

Interestingly, the branching ratios for iJ — > 65 and 
A ^ bb remain the dominant channels decay modes for 
modest (or larger) tan/3 even when decays to gauge bo- 
son pairs becomes kinematically accessible. For itla ^ 
mz, the mixing angle tan 2a — > tan 2/3, and thus H is 
mostly H^. Larger tan/3 implies {H^) <C (i?"), and 
thus the 3-point couplings HWf^W^^, HZ^Z^ are sup- 
pressed. Analogously, since there is no expectation value 
for the CP-odd scalar, these 3-point couplings are ex- 
actly zero. Thus, the decays into bb remain dominant 
until rriH.A ^ min(2/i, 2Mi, 2M2), where decays into the 
lightest gauginos becomes kinematically accessible. This 
suggests that the _ff , A — bb mode is viable up to well 
past 200 GcV (twice the smallest allowed Higgsino mass) , 
and for the Study Points in this paper, up to and beyond 
300 GeV. 
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FIG. 4. The branching ratio for decays to a Higgs boson 
is shown as a function of Mi/p for ~ 150 GeV, |^| = 
150 GeV, and tan /3 = 4. The upper plot shows the decay 
rates of heavy gauginos into the lightest Higgs boson, while 
the lower plot shows the summed decay rates to the heavier 
Higgs bosons H/A. The squark and slepton masses are taken 
to be 1 TeV. 



IV. MIXED HIGGSINO/BINO AS DARK 
MATTER 



One of the more attractive features of the weak scale 
supersymmetry with conserved _R-parity is that there ex- 
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FIG. 5. The branching ratio for squark decays to a Higgs bo- 
son as a function of Mi/fi for rriA = 150 GeV, \fi\ = 150 GeV, 
and tan /3 = 4. The upper plot shows the decay rate to the 
lightest Higgs boson, while the lower plot shows the summed 
decay rate to the heavier Higgs bosons H/A. As in Figure. [2] 
g_L refers to the sum of ul and d^, while qn refers to either ur 
or dfl.The squark and slepton masses are taken to be 1 TeV. 
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FIG. 6. The LSP reUc density as a function of h'h/^ with 
M2 = 2Mi, = 150 GeV, and tan/3 = 10. The squarks, 
sleptons, and tua were taken to be 1 TeV. The thermal relic 
density was calculated using micrOMEGAs v2.4 |33j . 



ists a stable, neutral, colorless, weakly-interacting par- 
ticle near the electroweak scale. In the post-LEP era, 
however, the prediction of present dark matter density 
does not automatically agree with the observation. LEP 
bounds typically forces the superpartner spectrum to 
be heavier and hierarchical. In this scenario, neutrali- 
nos are closer to pure gauge eigenstates, namely bino, 
wino and Higgsinos. Avoiding coannihilation and Higgs 
pole regions, the relic density is generally too large for a 
bino and too small for Higgsinos and winos. Matching 
cosmological data seemingly requires rather precise re- 
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lations among supersymmetry breaking parameters (see 
e.g. (34i). 

Much of these constraints follow from the requirement 
that the LSP relic density matches the observed cosmo- 
logical dark matter density. For our purposes, we are 
content to simply not predict too much dark matter, 
since other non-thermal sources of dark matter may be 
present. All of the Study Points considered in the paper 
automatically have a thermal LSP relic abundance that is 
at or smaller than the observed cosmological abundance, 
n^oh"^ < 0.1. 

In Fig. |6]we show the calculated thermal relic density 
U-^o is plotted as a function of Mi //i for a fixed value of 
1^1 = 150 GeV and tan /3 = 10. The thermal relic density 
was calculated using micrOMEGAs v2.4 [33]. All squarks 
and sleptons were taken to be 1 TeV. This clearly shows 
the variation of Xi relic density as the gaugino/Higgsino 
content of LSP is changed. For large values of Mi//i, the 
lightest ncutralino is mostly a Higgsino. As is well known, 
Higgsino-like neutralinos annihilate efficiently into gauge 
bosons, causing the calculated relic density to be smaller 
than the cosmological density. As the bino fraction of Xi 
increases with decreasing Mi, the annihilation rate goes 
down, and thus relic density goes up. Since the squarks 
and sleptons are much heavier than the gauginos, the 
bino rarely annihilates through them. For the specific 
parameters we considered, we find the annihilation rate 
can be optimized to give the right relic abundance to 
match the observed cosmological abundance when Afi ~ 

Matching the thermal relic density by taking Mi ^ 
means bino cannot decay into Higgsinos and Higgs 
bosons. Given t he n ear independence oi ql ^ W and 
qji ^ B (c.f. Sec III), only roughly half of the potential 
Higgs signal is lost given that right-handed squarks no 
longer lead to decays into Higgs bosons. We present one 
Study Points that demonstrates the Higgs signal remains 
perfectly viable when Mi ~ 



V. JET SUBSTRUCTURE ALGORITHM 

There are now several interesting techniques that ex- 
ploit jet substructure to enable better identification of 
standard model or beyond-the-SM signals P^SHTOl [531- 
1251 HZl [Ml [33 [35H37]. The central idea motivating the 
elaborate jet manipulation is that it is possible to seek 
a single "fat jet" (that is, a jet with a particular struc- 
ture consistent with one coming from a massive parti- 
cle decay) from the decay products of a boosted parti- 
cle. Butterworth, Davison, Rubin, and Salam (BDRS) 
^^ demonstrated that the Higgs boson of the Standard 
Model could be found with high significance using this 
technique [T3j . Their particular study has been validated 
by a realistic simulation done by the ATLAS collabora- 
tion |20!. 

The substructure algorithm developed by BDRS [13] 
to find a Higgs boson has two distinct parts: First deter- 



mine whether a jet contains substructure consistent with 
coming from a Higgs decay to bb. If it passes the criteria, 
"filter" the jet, improving the resolution of the invariant 
mass of the candidate resonance jet significantly. In order 
to identify a jet as a "fat jet" , BDRS stipulate two con- 
ditions: the mass of individual subjets are significantly 
smaller than the mass of the jet (the mass-drop condi- 
tion) and the splitting of the jet into the two subsets is 
not too asymmetric. The mass-drop condition basically 
checks how the jet-mass is distributed in the jet-area, 
and seeks out a jet that is consistent with one accom- 
modating all the decay products of a massive particle. 
Given the immense rate for QCD jets, the mass-drop 
condition alone is not enough. The background jets are, 
however, dominated by gluon splittings which exhibit soft 
and coUinear singularities. These singularities imply the 
majority of QCD subjets are asymmetric, so by rejecting 
particularly asymmetric splittings, the background can 
be further suppressed. 



Substructure for Supersymmetry 

In Ref. [H] we proposed an algorithm to extract 
a Higgs boson signal using its dominant decay mode, 
h ^ bb from a new physics event sample. Our algo- 
rithm exploits the techniques developed by BDRS, with 
some additional steps designed to allow our algorithm 
to be somewhat more efficient than BDRS when ap- 
plied to busier final states characteristic of new physics. 
Following the criteria laid out in Sec. |III[ the simplest 
superpartner cascade which yields a boosted Higgs is, 
q — )■ X + i — ^ x' + ^ + J: which necessarily involves one 
additional hard parton. More complicated signal events, 
with multiple extraneous, hard partons are easy to imag- 
ine. These hard partons, and their associated showers, 
can end up in the same fat-jet as the h — > bb. As these 
contaminating partons come from heavy particle decay 
and not from QCD radiation, they can survive the mass- 
drop and the asymmetry cuts (top and bottom quarks 
coming from the decay of superpartners are particularly 
dangerous as they also possess heavy flavor). Conse- 
quently, while declustering a fat jet, one may encounter 
multiple stages (say, "thresholds") that would pass all 
substructure criteria cited above. 

The BDRS algorithm is designed to consider only the 
first declustering stage that satisfies the mass-drop and 
asymmetry conditions, and as such, it is more suscepti- 
ble to false thresholds encountered in new physics events. 
BDRS jets are built using the inclusive C/A algorithm 
[35l441j . where subjets closer to each other are combined 
earlier, so the first threshold encountered will be where 
the subjets are maximally separated. To help distinguish 
between real and false thresholds, we need to use more 
information about the subjets. Although the contami- 
nating hard partons are not removed by the mass-drop 
and asymmetry cuts, they necessarily introduce a new 
scale into the jet. Rather than select a threshold based 
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on separation alone, we select the threshold where the 
subjet kinematics are maximally similar. Specifically, we 
impose a measure of similarity: maximize subjet hard- 
ness weighted by the inter-subjet separation. This mea- 
sure takes advantage of the isotropic decay scalar parti- 
cles - the Higgs bosons ~ that we are interested in. The 
algorithm [14J is described in full detail below. 

A. Our Algorithm 

The first step in our algorithm is to group final state 
particles, after all showering and hadronization, into 
"cells" of size Arj x Acj) = 0.1 x 0.1. All particles in a cell 
are combined, and the three-momentum of the total is 
rescaled such that each cell has zero invariant mass [24] . 
Cells with energy < 1 GeV are discarded, while the rest 
are clustered into jets. The initial clustering is done 
using the inclusive C/A algorithm, as implemented in 
FastJet |42], and taking the jet size to be i? = 1.2. Once 
the jets are formed, we search for heavy flavor; this is 
an essential step given that we want to discover the 
Higgs through its decays to bottom quarks. We b-tag 
jets by looking through the event record for b-mesons 
or b-baryons. If there is a b-flavored object within 20° 
of a jet direction, we tag the nearby jet as a b-jet with 
60% probability. If there are no b-flavored objects in 
the vicinity of the jet in question, the jet is tagged as 
a b-jet with a "fake-rate" of 2%. Every b-tagged jet in 
the event is then decomposed to search for substructure 
following the steps below: 

1. Undo the last stage of jet-clustering. As a jet is built 
from a sequence of 2 — > 1 mergings, unclustering one 
stage yields two subjets. The two subjets ji and j2 are 
labeled such that nij^ > rrij^ . 

2. Following Ref. [T3], subjets are checked for the ex- 
istence of a significant mass drop {rrij-^ < ii rrij ) , as 
well as non-existence of an asymmetry defined hy y = 

V l^m . > We use ^ ^ 0.68 and 

2/cut = (0-3)^ identical to Rcf. [13 . Both subjets are 
required to be 6-tagged and have pt > 30 GeV. If these 
conditions are satisfied, this stage of clustering (say, i-th) 
is recorded and then the following is calculated: 

min ) 

S. - , ^ i ^Rnn ■ (17) 

[PT,, +PtJ 

The quantity Si (namely, similarity) is an indicator of 
the similarity of the two subjets and is weighted by their 
separation ARj-^j^ . 

3. Replace j by ji and repeat from step 1 as long as j 
has further subjets. 

4. Select the stage of clustering for which Si is the 
largest. We anticipate that the two fe-tagged subjets, 
at this stage, are most likely to have originated from 
Higgs decay since they are more likely to be similar to 



each other. If the two C/A 6-tagged subjets originate 
from Higgs decay, the subjets with opening angle ARj^j^ 
should contain all the perturbative radiation from the 
hb system by virtue of angular ordering [43]. However, 
the subjets still tend to include too much contamination 
from underlying events. We then filter [13] the events: 
we cluster the jet constituents again using a finer angu- 
lar scale specific to the jet [we use, min (_/?;,(, /2, 0.3)] and 
retain only the three hardest components {bhg). Finally, 
we combine the three subjets and call the resultant a 
"candidate resonance jet" . 



B. Comparison with BDRS Algorithm 

Our algorithm declusters the event entirely, thereby 
checking multiple thresholds, while the BDRS algorithm 
only checks a single threshold. In an environment where 
there are few extraneous partons flying around, such as 
W^H production or even supersymmetric Higgs produc- 
tion from short cascades, there are few false thresholds 
and the two algorithms perform comparably. However, 
as the number of extra partons (and thus the number 
of false thresholds) increases, there is a clear difference 
in efficiency. Any threshold, genuine or not, will stop 
the BDRS algorithm, while our approach takes in all 
thresholds and sorts them out using the px similarity. 
Events with a true threshold masked by a false threshold 
at larger R will be missed by BDRS, but captured by 
our approach. Of course, the similarity variable will not 
always select out a true threshold from among several, 
so accuracy is not necessarily increased. Fig. [7] shown 
below, is a simple demonstration of how our algorithm is 
more efficient in a crowded environment. 

This Figure was generated from a signal sample of 80K 
PYTHIA-generated events using the spectrum SHSP 3 
listed in Table |TJ The key feature of this spectrum is that 
the gluino is as light as the squarks, so the signal sam- 
ple has a significant number squark-gluino associated- 
production events. The gluinos decay through off-shell 
squarks, and typically lead to busy, multi-jet events. As 
a measure of the increased efficiency, we can count the 
number of event under the putative Higgs peak and com- 
pare with the number of events in the bins adjacent to 
the peaks for each of the algorithms. The significance, 
defined as (# events in peak above continuum)/^(# con- 
tinuum events) is larger for the similarity algorithm by 
a factor of ~ 1.2. As we will later see, the continuum 
supersymmetric events are often the largest background 
to the Higgs peak, so we can expect the complete signifi- 
cance (including SM backgrounds) to increase by roughly 
the same amount. 

It is important to point out that, although we employ 
this algorithm to find Higgs, all we really check for is 
a massive particle decaying to 2 6 partons. Any heavy 
multiplet that decays to 2 6 should also be selected by 
our jet algorithm as long as they are boosted. Among 
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A. Supersymmetric Higgs Study Points 

The efficiency of our algorithm to find Higgs bosons is 
demonstrated on a set of benchmark points, Supersym- 
metric Higgs Study Points (SHSPs), tabulated in Table|l] 
These Study Points are grouped into three categories. 

• Study Points 1,2 represent spectra in the decou- 
pling limit {mA = 1 TeV) . In SHSP la and lb the 
LSP is mostly bino, all squarks are at 1 TeV, and 
the sleptons are at 1 TeV and 350 GeV respectively. 
In SHSP 2a and 2b the LSP is a maximal mixture of 
Higgsinos and bino. In SHSP 2a once again we use 
heavier squarks and sleptons while slightly lighter 
squarks and sleptons are used in SHSP 2b. 



FIG. 7. Comparison of the candidate resonance jet mass us- 
ing the similarity algorithm (black) and the original BDRS 
algorithm (red), applied to our ^-tagged fat jet sample. The 
signal point used for the comparison is SHSP 3 (see Table |l|, 
and the vertical axis has been rescaled to correspond to an 
integrated luminosity of £ = 10 fb~^. While the accuracy of 
the two algorithms is similar, the similarity algorithm is more 
efficient. 



SM particles we expect to find Z ^ bb^ Also, when 
both H and A are light and decay to bb, our algorithm 
can discover Higgs bosons as long as they are produced 
in a superpartner cascade. 



VI. RESULTS 

Having demonstrated sparticle cascade decays as a vi- 
able, important source of boosted Higgs bosons and de- 
scribed our substructure algorithm in detail, we now 
demonstrate the effectiveness of our proposal. To best 
convey our results, we first propose a collection of Study 
Points on which we use the candidate resonance jet find- 
ing algorithm. While by no means exhaustive, the Study 
Points have a diverse set of MSSM parameters. After in- 
troducing the Study Points, we then list the set of back- 
grounds we considered for this work and show the way 
in which sets of conventional cuts can be used to reduce 
these. The candidate resonance jet finding algorithm is 
then run on this set of rarefied events (both signal and 
background events). Finally, masses of the candidate res- 
onance jets are plotted to estimate the signal significance. 



^ In practice, as light jets can occasionally fake b-jets, any boosted, 
heavy particle which decays hadronically (t, W, ■ ■ ■ ) has a chance 
of being picked up by the substructure algorithm. 



• SHSP 3 has Mi ~ and large iriA, such that 
the LSP has a thermal relic density that matches 
cosmological measurements. 

• The final set of Study Points, SHSP 4,5 = 
150 GeV) and SHSP 6 (|/i| = 200 GeV) are repre- 
sentatives of spectra in the smaller niA region. The 
main difference between SHSP 4,6 versus SHSP 5 is 
the sign of the term. As shown in Fig. [5] when 

is low the sign of /i greatly influences which Higgs 
bosons the gauginos decay into. For SHSP 4, 6, de- 
cays to h predominate, while H/A predominate in 
SHSP 5. 

To simulate the supersymmetric signal, we use 
PYTHIA v6.4 to generate parton level events, with sub- 
sequent showering and hadronization. The lowest-order, 
inclusive superpartner production cross sections are large 
(©(pb)) and are listed for all Study Points in Table |T] 
These cross sections are somewhat misleading, since the 
quoted cross sections also include electroweak production 
of light charginos and neutralinos. In the scenarios we are 
considering, the lightest charginos and neutralinos have 
a large Higgsino component and thus large couplings to 
the Z boson. As a result, the LHC cross sections for 
neutralino pair production XiX2j chargino pair produc- 
tion xfxf Eind associated chargino-neutralino produc- 
tion X? 2X1^ cire all quite large, 0(0.5 — Ipb). While a 
large chargino/neutralino production cross section will 
likely enable the discovery of new physics, light neu- 
tralinos and charginos do not decay to Higgs bosons so 
these events are of no use for a Higgs search. There- 
fore, in order to fairly judge our Higgs-finding algorithm, 
we have included the fraction of supersymmetric events 
containing a Higgs boson (h/H/A) in Table |lj This frac- 
tion was calculated by counting the number of on-shell 
Higgs bosons, without any kinematic cuts, in samples 
of PYTHIA-generated supersymmetric events. The final 
row in Table [ij (Jh/H/A, is simply the inclusive supersym- 
metric cross section times the fraction of supersymmetric 
events containing a Higgs boson. 
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SHSP la / SHSP lb 


SHSP 2a / SHSP 2b 


SHSP 3 


SHSP 4 


SHSP 5 


SHSP 6 


tan /J 


10 


10 


10 


5 


6.5 


10 


Ml 


300 GeV 


150 GeV 


163 GeV 


200 GeV 


200 GeV 


300 GeV 


Ma 


600 GeV 


300 GeV 


400 GeV 


400 GeV 


400 GeV 


600 GeV 


Ms 


2.1 TeV 


1.05 TeV 


1.0 TeV 


1.4 TeV 


1.4 TeV 


2.1 TeV 


fJ- 


150 GeV 


150 GeV 


200 GeV 


200 GeV 


-150 GeV 


150 GeV 


niA 


1 TeV 


1 TeV 


1 TeV 


150 GeV 


150 GeV 


200 GeV 


at 


900 GeV 


-900 GeV 


900 GeV 


2.04 Te\f] 
1 TeV 


1.4 TeV 


900 GeV 


ruq 


1 TeV 


1 TeV/750 GeV 


1 TeV 


1 TeV 


1 TeV 


mj- 


1 TeV/350 GeV 


1 TeV/350 GeV 


350 GeV 


1 TeV 


1 TeV 


1 TeV 


mh 


116 GeV 


117 GeV 


116 GeV 


114 GeV 


115 GeV 


115 GeV 


run 


1 TeV 


1 TcV 


1 TeV 


161 GeV 


157 GeV 


202 GeV 


niA 


1 TeV 


1 TeV 


1 TeV 


150 GeV 


150 GeV 


200 GeV 


mH± 


1 TeV 


1 TeV 


1 TeV 


169 GeV 


170 GeV 


216 GeV 


Xi 


138 GeV 


110 GeV 


140 GeV 


157 GeV 


136 GeV 


138 GeV 


X2 


-158 GeV 


-161 GeV 


209 GeV 


-207 GeV 


-163 GeV 


-158 GeV 


X3 


206 GeV 


174 GeV 


-209 GeV 


227 GeV 


210 GeV 


306 GeV 


X4 


625 GeV 


338 GeV 


429 GeV 


433 GeV 


426 GeV 


623 GeV 


+ 

Xi 


148 GeV 


137 GeV 


191 GeV 


187 GeV 


152 GeV 


148 GeV 


x^ 


625 GeV 


337 GeV 


429 GeV 


433 GeV 


426 GeV 


623 GeV 


O"tot 


3.9 pb 


5.8 pb / 8.07 pb 


2.76 pb 


2.4 pb 


4.1 pb 


4.0 pb 


% Higgs 


4.5%/3.4% 


4.2%/6.8% 


6.6% 


12.8% 


8.6% 


7.0% 


O'h/H/A 


0.18 pb/0.13 pb 


0.24 pb/0.55 pb 


0.18 pb 


0.31 pb 


0.35 pb 


0.28 pb 



^ = at for this point as well 

TABLE I. The parameters, part of the spectrum, and some relevant coUider information for the Study Points used in this 
analysis. The spectrum was computed with SUSPECT2 |44) . The quoted cross section is determined at lowest order for the 
LHC operating at a center of mass energy of ^/s = 14 TeV. See the text for the definition of % Higgs and (t^iuia- 



Backgrounds and Cuts: 

The primary SM backgrounds we consider are: 

• tt+jets 

• heavy flavor 

• WjZ^ jets 

• tt + hh 

These backgrounds are familiar from many supersym- 
metry/BSM searches. They have large cross sections, 
multiple jets, some of which are 6-jets, and sources of 
missing energy from vector boson decays. The back- 
ground events are first generated at parton-level using 
ALPGEN vl3 [45] and are then showered and hadronized 
using PYTHIA v6.4[^ We also use the ATLAS tune [15] 
in PYTHIA to model the underlying event. Jet manip- 
ulation is done using Fast Jet [32]. We do not perform 
any detector simulation or smearing of jets. A realistic 
ATLAS/CMS specific search in the spirit of Ref. 20, is 
beyond the scope of this work. However, since high 
jets result in a large amount of energy deposited in the 
calorimeter cells where energy resolution is excellent, we 



3 All events generated with ALPGEN using CTEQ5L par- 
ton distribution functions and default options for factoriza- 
tion/renormalization scheme. 



do not expect smearing to significantly modify our re- 
sults. 

Before we run our substructure algorithm, we intro- 
duce cuts to isolate the signal from the background. 
Rather than tailoring the cuts to each specific SHSP 
point, we choose a more generic set which can be applied 
to all Study Points. In particular, we use: 

\. $T> 300 GeV. 

2. S"*" jets, at least one of which is tagged as a b-jet. 
To be counted as a jet, we require pT > 200 GeV 
- the "Pt requirement on the jets is set so high be- 
cause we want to capture an entire boosted object 
(ideally a Higgs) within a single jet. As explained 
in Sec. |Vj objects will be reconstructed from sub- 
jets contained within individual high-p^ jets rather 
than combining multiple jets. We impose a pseu- 
dorapidity cutoff of \r\\ < 4.0 for jets which are not 
flavor-tagged, while b-tagged jets are restricted by 
the pseudorapidity extent of the tracker, |?7| < 2.5. 

3. No isolated leptons with pr > 20, |?7£| < 2.5. 

4. ijj^'" = J2,PTA > 1-0 TeV, where the sum extends 
over all jets indexed by i. 

Large missing energy, large Hi^^^, and high jet multi- 
plicity are often the characteristics of new physics and, 
in particular, of weak scale supersymmetry with R- 
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parity |^ These variables are widely used in supersym- 
metric searches and we use them here. After I^t and 
H^^^^ cuts, the biggest background is by far tt+jets. In 
order to suppress the tt+jets further we introduce a lep- 
ton veto; the logic behind this cut is that any tt+jets 
events which pass the large I/^t cut most likely contain 
at least one leptonic . 

We collect all events that pass our preliminary cuts 
and run the substructure algorithm described in Sec. |V] 
Events which pass the substructure selection have at least 
one 6-jet with substructure and, consequently, at least 
one candidate resonance jet. 

The assumed background cross sections and their efh- 
ciencies under the imposed cuts are summarized below in 
Table [n] To show how substructure cuts affect the signal 
and background, we have broken up the efficiencies into 
two stages. The first stage, Ccuts, is calculated after the 
'conventional' cuts - ^T,Hi^^^, jet multiplicity and lep- 
ton veto - have been imposed. Then, after running the 
substructure algorithm, the surviving events are counted 
to determine ecuts+subs- 



Process 


(7 


^cuts 


^cuts-\-suhs 


tt+ jet 
tt+ 1 jet 
tt + 2+ jet 


474 pb 
248 pb 
132 pb 




9.2 X lO"'' 
2.1 X 10"* 




1.62 X 10"" 
4.84 X lO"'^ 


tt + bb 


0.83 pb 


1.9 X 10"* 


4.6 X 10"'' 


W{£u) + 2 jets 
W{eiy) + 3+ jets 


127 pb 
50 pb 


2.3 X 10"" 
2.3 X 10~* 




1.08 X 10"'^ 


Z{w) + 2 jets 
Z{w) + 3+ jets 


80 pb 
29 pb 




2.3 X 10"* 




1.47 X 10"'' 


Z{uiy) + bb 
Z{vu) + &6+jet 


1.4 pb 
1.4 pb 




3.5 X 10"* 




6.94 X lO"'' 


W{lv) + bb 
W{lv) + bb+iet 


1.1 pb 

2.2 pb 


2.6 X 10"'' 
1.1 X 10"* 




3.6 X 10"'^ 


SHSP lA 


3.92 pb 


0.015 


3.7 X lO"'^ 



TABLE II. Signal lA and background cross sections and effi- 
ciencies. The efficiency for the other Study Points is similar. 
All backgrounds were generated with parton level cuts on jets: 
PT,j > 30 GeV (25 GeV for tt + jets), \rij\ < 4.0, ARjj > 0.4. 
An additional cut > 75 GeV was used for all W/Z back- 
grounds, while \rit\ < 2.75 was added for all heavy- flavor 
events. All background cross sections are LO except for 
tt + jets, which has been rescaled to the NLO MCFM gT] 
result a = 855 pb {K ~ 1.8). 

As can be seen from the Table, the conventional cuts 
are quite effective at reducing the background. The 
signal efficiency under the conventional cuts looks low. 



However, as explained in Sec. VIA many supersymmet- 
ric events for these Study Points come from electroweak 



Due to our large jet-py requirement, H!^'^^ calculated with our 
definition can be quite different than the sum of all visible trans- 
verse energy in the calorimeters (often referred to as Hx cal)- 
However, Hj- ^^l > which we would rely on for triggering, will al- 



ways be bigger than H; 



■jets 



chargino/neutralino pair production which do not con- 
tain the sufficient energy or jet multiplicity to pass our 
cuts; the efficiency for squark/gluino initiated events is 
higher. Requiring jet substructure suppresses the back- 
ground further relative to the signal, however the real 
power from substructure comes in the shape of the jet- 
mass distribution. Therefore, the final step in our search 
strategy is to plot the invariant mass of all candidate res- 
onance jets and look for a peak consistent with a Higgs 
boson. The candidate resonance jet mass plots for each 
of the 8 benchmark Study Points in Table |T] are presented 
in the following sub-sections. To break up the results, we 
have grouped the Study Points into the same three cate- 
gories used in Table |l] high-mA, low-m^, and one Study 
Point with a LSP thermal relic density that matches cos- 
mological observations. 

In all of the following plots, the contribution from all 
supersymmetric events (inclusive superpartner produc- 
tion) are shown together on top of the SM background. 
While the supersymmetric contribution contains our sig- 
nal, Higgs bosons from sparticle decays, it also contains 
new backgrounds. Top quarks and W/Z bosons will also 
be copiously produced in cascade decays and can occa- 
sionally pass the substructure cuts. In fact, in several 
circumstances this supersymmetric background is larger 
than the SM backgrounds. 



B. High-rriA points: SHSP la/lb, 2a/2b 

The first of the high-m^ points, SHSP la and 16, are 
characterized by small /i. The large mA kinematically 
forbids squark decays to other Higgs states (H/A/H^), 
while the low /i = 0.5Mi — 150 GeV implies a very 
Higgsino-like LSP and thus large branching fractions 
— > /i + xfiXsA h + Xi.2- These points are ide- 
ally suited to our analysis, and the resulting candidate 
resonance jet mass plot. Fig. [8] verifies this. The peak 
arising from Higgs decay is unmistakable over the rela- 
tively featureless SM background. 

The small shoulder to the left of the Higgs peak comes 
from Z bb events. Higgs bosons and the Z are pro- 
duced in relatively equal amounts for these two points, 
due to the arguments presented in Sec. |llj However, the 
Z -^bb branching fraction is only 1 /6 as large as /i — > bb, 
and the resulting Z peak is small. 

The only difference between SHSP la and lb is the 
mass of the sleptons. In SHSP lb, the sleptons are 
light enough that the heavier charginos and neutralinos 
can decay into them. As demonstrated in Figs. |1|2[ 
new chargino/neutralino decay modes imply a smaller 
fraction of decays to Higgs bosons, and thus a smaller 
signal. However, comparing the top and bottom 
plots in Fig. [Sj we can see that the rate decrease to due 
decays to sleptons of mass Mi < mj < M2 is quite minor. 

To get a quantitative idea of how well our algorithm 
can find the Higgs, we estimate the significance of the 
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FIG. 8. Distribution of the candidate resonance jet mass nor- 
malized to 10 fb~^ of integrated luminosity at 14 TeV center of 
mass energy. The contribution from supersymmetric particles 
is shown in blue for points SHSP la (top) and SHSP 16 (bot- 
tom). Standard Model contributions, which come primarily 
from tt-|-jets, are indicated by the red and green histograms. 



Higgs peak on top of the SM and continuum new physics 
background. We determine the SM plus continuum con- 
tribution using the same simple method as in Ref. '14]; 
the histograms 1 — 2 bins on either side of the Higgs 
peak are connected with a Hne, and anything within the 
resulting trapezoid is counted as background. To check 
the veracity of this procedure, we have looked back into 
the signal events and assigned each event with a candi- 
date resonance jet to an initial parent parton {t/W/Z/h) 
according to which heavy particle is closest in AR We 
find the fat -jets with a parent Higgs are indeed confined 
to the peak and neighboring bins. Events with a Z par- 
ent are similarly confined to the bins near mz, while the 
continuum events are composed of t/W events. Using 
the peak — 1/ + bins to define the signal region (mean- 
ing the bins —2/ + 1 on either side of the peak are used 
to determine the background), we find S = S/Vb = 7.9 
for point SHSP la. The same procedure, taking the sig- 
nal region to be the peak ±1 bins, gives S = 9.6 for 
point SHSP 16. These significances are just rough esti- 
mates. We have taken quite aggressive conventional cuts 
to render the SM background as small and featureless as 
possible; less strict cuts may lead to higher significances. 



as would optimization of the cuts for each SHSP point. 

The next two high-my^ points are more challenging for 
three reasons. First, points 2a and 26 have a smaller 
Ml. As we saw in Figs. T]2 a lower Mi/ jj. means fewer 
Higgs bosons are produced from squark cascades. Sec- 
ond, lowering Mi while holding Mi : M2 : M3 ratio 
fixed implies a much lighter gluino. While the gluinos 
are light in this scenario, they are still capable of decay- 
ing to on-shell squarks, so supersymmetric events orig- 
inating from gluinos - either from gluino pair produc- 
tion or squark-gluino associated production - have more 
jets than events originating from squark pairs. Addition- 
ally, because gluinos decay democratically to all species of 
squarks, gluino cascades can easily include top and bot- 
tom quarks. The third difficulty with 2a and 26 comes 
from right-handed squarks. Right-handed squarks, pro- 
duced either in pairs or associated with a left-handed 
squark or gluino, decay to bino plus jet, with the bino 
in this spectra spread between Xi,X2, and X3. However, 
as can be seen from Table. [H] the mass-gaps among the 
three lightest neutralinos are small enough so that most 
two-body decay modes are kinematically forbidden; the 
neutralinos decay instead via an off-shell W/Z/h plus a 
lighter chargino/neutralino. Off-shell, hadronic decays 
quickly lead to an increase in the number of hard par- 
tons in the event. For example, a typical signal process 
involving one right-handed squark: pp — > qhqR followed 
by qh xl+i Xi + h + j and qji x" + j ^ X? + 3j 
involves 4 extra hard partons, any one of which can fall 
in the same fat-jet as the Higgs boson. Longer cascades, 
coming from gluino production or more decay steps, are 
easy to imagine and will contain even more hard partons. 

When extra partons from superpartner decays are erro- 
neously combined with all or part of a Higgs candidate, 
the jet mass becomes smeared. The smearing is exac- 
erbated by the fact that, following BDRS, we take the 
three hardest subjets during filtering to define the can- 
didate resonance jet. Such a procedure remarkably im- 
proves the mass resolution of a Higgs jet when the correct 
threshold is identified and none of the extra hard partons 
produced in association with the Higgs is inside the Higgs 
cone. The three hard partons during filtering then cor- 
rectly capture 66 from Higgs as well as the first radiation 
inside the 66 system. On the other hand, if there is an 
extra hard parton inside the 66 system, the filtered reso- 
nance jet may end up containing 66 + hard parton instead 
of being 66 + radiation and consequently having a skewed 
invariant mass. This smearing is clearly visible in Fig. [9] 
and creates the feature extending from the Higgs peak 
to higher mass. However, despite all the contamination 
from auxiliary hard partons, the Higgs peak is still quite 
visible. Perhaps more elaborate subjet algorithms could 
be used to clean up the high-mass tail further. 

Moving from 2a to 26, the squark mass decreases. 
Lighter squarks are produced even more prodigiously, 
as reflected in the enormous superpartner cross section, 
however they impart a smaller boost on their decay prod- 
ucts. The increased production of right-handed squarks 
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FIG. 9. Distribution of the candidate resonance jet mass in 
points SHSP 2a (top) and SHSP 2b (bottom). As in Fig. [s] we 
assume 10 fb^^ of integrated luminosity and a 14 TeV center 
of mass energy. 



in SHSP 2b, a factor of ~ 4 compared to point 2a is 
responsible for the increased number of supersymmetric 
events away from the Higgs peak. The slepton mass in 2b 
is also smaller than in 2a, however this has only a small 
effect since the sleptons are still too heavy for the higher- 
tier charginos and neutralinos to decay into, mj > M2- 
Estimating significance in the same way as we did for 
SHSP la/6 and using — 1/ + 2 bins to define the signal 
region, we find a significance of (3.8, 5.6) for points SHSP 
(2a, 26). 

Having seen the effects of decreasing the squark mass, 
it is natural to ask what happens if we do the opposite 
and raise niq and M3 while keeping the rest of the super- 
symmetry parameters fixed. The squark/gluino mass sets 
the scale for the boost of its subsequent decay products, 
including any Higgs bosons. One may worry that a higher 
sparticle scale would lead to Higgs decay products which 
are so boosted that finite detector granularity or the need 
to b-tag multiple subjets would render our algorithm use- 
less. This does not occur, however, as is evident in the 
distribution of the subjet angular scale Rbb- We find a 
rather flat distribution between 0.3 < Rbb < 1-2, which 
persists even as mq/M^ is raised to several TeV (squark 
with mass beyond 3 TeV have such low production cross 



section that they become phenomenologically irrelevant 
at the LHC). Therefore, even if granularity/tagging inef- 
ficiencies ruin the most highly boosted Higgs bosons, the 
broad tail of Rht indicates that our algorithm can remain 
viable throughout the range of interesting squark masses. 



C. Relic Point: SHSP 3 

The parameters of point SHSP 3 have been chosen such 
that the LSP has a thermal relic abundance that matches 
cosmological observations for the dark matter density. 
As described in Sec. |IV[ this requires delicately adjust- 
ing Ml ~ \fj,\ to get the right bino/Higgsino admixture 
in the LSP. Point SHSP 3 also has light sleptons, how- 
ever the branching fraction to Higgs bosons is still high 
enough for our analysis to be successful. The cascade of 
superpartner decays contain more Z than h, due to the 
smaller value of Mi , but the Higgs peak remains clearly 
visible. 




80 100 120 140 160 180 200 220 240 

candidate resonance jet mass (GeV) 

FIG. 10. Distribution of the candidate resonance jet mass in 
point SHSP 3. The parameters for this point were chosen to 
produce the correct dark matter relic abundance for the LSP. 
As in Fig. [s] we assume 10 fb"'^ of integrated luminosity and 
a 14 TeV center of mass energy. 



D. Low-mA points: SHSP 4-6 

The final set of Study Points have smaller ttia and a 
small value for tan/3. The region of smaller uia, small 
tan /? is known to be difficult for traditional MSSM Higgs 
boson searches, so these points serve as an important test 
of our algorithm. To ensure the lightest Higgs boson has a 
mass that exceeds the LEP bound, we allow larger mixing 
in the stop (and sbottom) sectors. These Study Points 
are therefore quite similar to the "maximal-mixing" sce- 
nario often considered in collider searches j48j . 

Perhaps the most interesting consequence of tua <^ 
mg,M2, is that the heavier Higgs bosons H/A also ap- 
pear in the superpartner decay cascades. For mHymA — 
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150 GeV, the H/A decay predominantly into bh and are 
light enough that they will emerge from sparticle decays 
carrying a substantial boost. With these characteristics, 
H/A will be captured by our algorithm. This opens the 
exciting possibility, shown in fig. (VID), of discovering 



multiple distinct Higgs bosons with a single analysis. 
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captionDistribution of the candidate resonance jet mass in 
points SHSP 4 (top), SHSP 5 (middle), SHSP 6 (bottom). 
As in Fig. ([sjl we assume 10 fb~^ of integrated luminosity 
and a 14 TeV center of mass energy. 



In SHSP 4, the top plot of Fig. |VI D| heavier charginos 
and neutralinos decay to h rather than H/A making the 
h peak unmistakable. Some H/A are present, and lead 
to the feature near tua = 150 GeV. Given the size of the 
iriA feature and its proximity to the top mass, detector 
resolution effects, which we have treated very simply in 
this paper, become more important and need to be taken 



into account correctly. H/A discovery will likely require a 
more specialized analysis, but it is certainly possible that 
both H/A and h could be discovered with this technique 
given sufficient integrated luminosity. 

In point SHSP 5, the /i term is negative. With /i 
and mA similar in magnitude, the Higgs mixing ma- 
trix becomes particularly sensitive to the relative sign 
between these two mass parameters and cancellations 
can occur once couplings are expressed in terms of mass 
eigenstates. For /i < 0, the h coupling to higher-tier 
charginos/neutralinos is suppressed by one such cancella- 
tion, and cascade decays to H/A are more likely. We can 
clearly see this effect in fig. ( | VI D I ; the h peak is barely 
visible over the continuum new physics events, while the 
narrow H/A peak at 150 GeV is clearly evident. 

The final point, SHSP 6, has exactly the same su- 
persymmetry parameters as SHSP la except niA — 
200 GeV. This is the ideal point for detecting both 
the light and heavy Higgs bosons with a single analysis. 
The niA is low enough that Xi X2 have a moderate 
branching ratio to H/A, while ttia is heavy enough to 
avoid getting mistaken for new physics continuum or a 
top quark. Taking the signal region to be — 0/ -f 1 bins 
(—1/ -|- 1) around the h peak, we find a significance of 
(3.9,8.2) for points (SHSP 4, SHSP 6). Repeating the 
same procedure around the H/A peak, we find a signif- 
icance of (5.2,4.5) for (SHSP 5, SHSP 6) using signal 
regions — 1/ + bin. 

Low values for imply light charged Higgs bosons, 
which are constrained by the fiavor process 6 — >■ s + j. 
While the specific spectra we are looking at have b — > 
s + "f slightly larger than the experimentally allowed 
range [331 135] , slight changes in the spectrum, such as 
lowering the third generation squark masses or introduc- 
ing squark mixing can introduce cancellations and signif- 
icantly alter the branching ratio 5 — > s + 7 [50] . These 
changes to the spectrum need not effect the supersym- 
metric Higgs signal. Therefore, in the same spirit as [31I3], 
we focus on direct Higgs detection prospects and ignore 
indirect constraints for the time being. 



VII. DISCUSSION 

The power of using jet substructure with boosted Higgs 
decays into bb suggests the search for the MSSM Higgs 
bosons should be entirely rethought and redone, with full 
detector simulations. Our estimates, without jet energy 
smearing and without a realistic detector simulation, sug- 
gest that with less than 10 fb~^ of data at ^/s = 14 TeV, 
signal significance can exceed 5 for the h — > bb channel 
alone given total superpartner production rate of order a 
few pb. This is possible given the outstanding mass res- 
olution of our reconstruction technique combined with 
the power that jet substructure provides in discriminat- 
ing Standard Model and supersymmetric backgrounds. 
We have been relatively conservative in our candidate 
resonance jet finding algorithm given our fiat 6-tagging 
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efBciency: we required a triple b-tag - the original jet as 
well as two subjets. Nevertheless, our estimates of signal 
significance are just that - estimates. We urge the AT- 
LAS and CMS collaborations to carry out full detector 
simulations, along the lines of what was done by ATLAS 
to study the boosted Higgs into 66 mode in the Standard 
Model [20] ■ 

The notion that both h as well as H and A could be 
found using jet substructure techniques is particularly 
interesting given the difficulty that conventional search 
strategies have within the smaller niA and smaller tan/3 
region. The ATLAS and CMS TDR suggest fully cov- 
ering the MSSM parameter space requires considerable 
integrated luminosity, 60-100 fb~^. Our technique has 
the potential to cover this region much more rapidly. 

It is interesting that the MSSM parameter region most 
favorable to finding a signal of Higgs bosons is also the 
one with the least fine-tuning, namely, small /i (e.g. jSlj). 
Nevertheless, gaugino mass unification, and other as- 
pects of the superpartner hierarchy are somewhat less 
constrained. 

Finally, finding evidence for Higgs bosons within a new 
physics event sample provides an incredibly important 



connection between the new physics and the Higgs sector 
- i.e., the Higgs sector is necessarily coupled with the new 
physics. This connection can be established far faster 
than sorting out which kind of new physics is present 
based on the population of different BSM search chan- 
nels. The generic search strategy proposed here builds 
on our previous paper |14| . demonstrating the power of 
this method applied to the MSSM with a neutralino, or 
neutralino-equivalent, lightest supersymmetric particle. 
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